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Abstract We obtain some solutions to the Einstein-Klein- Gordon equations without a cosmo- 
logical constant for an exponential potential scalar field in a Bianchi VIo metric and investigate 
their behaviour. 



o. 

; I. INTRODUCTION 

OX) 

Inflationary theories are expected to explain the observed isotropy of the Universe by assuming an inflationary 
expansion in the early universe, since as the so-called "Cosmic no-hair" theorems state, all the solutions of the 
ly-j . Einstein equations initially expanding and with a positive cosmological constant evolve towards the de Sitter solution 
[l|, @. However, only those Bianchi models that have the FRW models as particular solutions isotropize Q. One of 
the potentials which has received special attention is that of Liouville form (exponential potential), the one we will 
'■ use in this work. 

Some anisotropic cosmological models have already been studied. Aguirregabiria, Feinstein and Ibanez Q analyzed 
the Bianchi I models with exponential potential by reducing the problem of finding exact solutions to the solution 
of one third order differential equation. In Q a non-local transformation was used to linearize this equation and the 
— —i- general solution was found. In this work we shall extend this analysis to the Bianchi VIo, reducing, again, the problem 
of finding exact solutions to resolve a nonlinear differential equation and investigating the behaviour of their solutions 
^ . at early time and its asymptotic stability in the far future. 

■ The line-element for a Bianchi VIo cosmological model can be written in the following form: 

§J ! ds 2 = e /(t) (-dt 2 + dz 2 ) + G(t) (e z dx 2 + e- z dy 2 ) . (I) 



and the corresponding Klein-Gordon and Einstein field equations for the metric (|TJ) are as follows: 



00 

*+%i + e?%r=0, (2) 



G 



% = 2e*V, (3) 



We will consider a homogeneous self-interacting scalar field (<f> = <fi{t)) with an exponential potential V — Ae k ^, 
which is the only potential that separates the Klein-Gordon equation in two parts, one containing only the scalar 
field and the other containing only geometrical quantities (see [a]). Making use of this potential we can express the 
scalar field as follows 



m k G 
G ~ 2G' 



(5) 
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After some algebra, the above equations can be reduced to the following one: 

G 2 G-KGG 2 -GGG+^GG 2 +m 2 G = 0: (6) 

where K = fc 2 /4 — 1/2. This is the equation we are going to study in detail. First we will consider the behaviour of 
the solutions when t — > and then we will study their stability. 



II. BEHAVIOUR OF THE SOLUTIONS WHEN t -> 



A. Particular case 



If we consider Equation ^ when m = 0, and we assume that the function G has the form G(t) = Go{At) n when 
At — > (At = t — to), we find two types of solutions, one with n = 1, which leads to a vanishing potential, and the 
other with n = 1/K. In the special case of — A < K < we obtain a power law inflationary solution (see Q), which 
can be written in its synchronous form as follows: 



ds 2 



dT 2 + PT^i e z/ ^dX 2 + e- z/ ^dY 2 + dZ 2 



(7) 



where p=[(n+ 2)/{2y^)} 2/{n+2) , a = and O a constant. 

We have studied the asymptotical behaviour of the scalar curvature, whose expression for the metric ^ is: 



R = e-W 



1 + 

2 2 G 2 G J 



(8) 



In this particular case, we obtain that R diverges at t — for K > 0, so we can speak about a singularity at that 
point. In the case of negative K the scalar curvature vanishes when t — > 0, due to the fact that t — > means T — > oo 
in Equation (J7]). 



B. General Case 



In the most general case the equation we have to deal with is: 

^ [„. - (1 + A>» + - (I + ^) ^ [„* - „] ^0. (9) 

When n < 1 we obtain the same results that we had in the particular case. Solutions with n > 1, however, are not 
permitted. This means that solutions with > K > 1 do not appear in the general case. 

As we did for the general case, we can calculate the behaviour of the scalar curvature, given by equation (j8]). The 
relevant solutions (those with — ^ <if<0orl< K) behave as they did in the m = case, and therefore, the 
conclusions about the divergence of R remain unchanged. 



III. STABILITY OF THE ASYMPTOTIC SOLUTION 



A. Particular Case 



To investigate the behaviour and stability of the solutions, we can rewrite Equation © in a quite different way, by 
introducing a new function h — ^ and redefining the time variable as follows: 

dr] = hdt = d\nG. (10) 

After some transformations Equation © reads 



d_ 

dr) 



h 12 h 2 1 

h K In ft, 

2 2 2 



+ <! + *) 



(11) 
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where primes are derivatives with respect to the new variable r\. This is the equation of motion for a dissipative or 
antidissipative system, with the potential V(ft) = — | In ft. 

Equation (fTTj) presents local minima when ft 2 , = , for K > (i.e. k 2 > 2). As the dissipative term given 
by the right-hand side of Equation (fTTj) is negative definite in the asymptotic regime, the bracket on the l.h.s. of 
Equation ([lTj) define a Liapunov Function Q , [i| , [l(| ■ Then, the corresponding exact solution, expressed by 

G = Goe^' (12) 

(with Go a constant) is stable for t — > oo and for any initial condition. This result allows us to study the behaviour 
of the solutions around these equilibrium points. To first order in perturbations and for K > — i, there is a two- 
parameter family of stable solutions that behaves as Equation (|12[) . The trajectories in the phase plane (ft, ft) can be 
divided in two different groups: for K > i the solutions cut the axis in the phase plane an infinite number of times, 
so they spiral around the constant solution fto (Equation p^|) ). For — | < K < g the solutions do not cut the ft axis 
or they cut it once. 



B. General Case 



Equation © when m = can be written using ft and rj variables, as follows: 



d_ 

dr) 



V(ft) 



= 2m 2 e- 2 "lnft- \ (1 + K)- 



1 + 2m 2 e- 2l) 
2ft 1 



-2.7 



} In ft, and the local minima will be given by ftp = - 



(13) 



+2m 2 e~ 2 " 



2K 



where now the "potential" is V(ft) = ^ - {5 + r 

One has to be more careful in this case because the equilibrium points are stable only if K > i as a consequence of 
Liapunov's theorem. In this case, the equilibrium points are not fixed, they depend on 77, and the equation that rules 
their behaviour gives us the following solution: 



G, ; 



± 



m 



(t-to) 



(14) 



so the final asymptotically stable solution behaves again as Equation (|12[) . which represents an anisotropic solution. 
This result had been obtained numerically in (Tl| and, in our case, it corresponds to the simplest solution we can 
obtain from Equation ([5]) when m = 0. 



IV. CONCLUSIONS 



We have studied the general behaviour of Bianchi VIo solutions with a self- interacting scalar field and an exponential 
potential, by reducing the system of field equations to one differential equation of third order, Equation ©. 

We have analyzed the behaviour of these solutions when t — > and have found several anisotropic solutions which 
appear to have a singularity at t = 0, since the scalar curvature diverges there. A power-law inflationary solution 
corresponding to fc 2 < 2 has also been found. 

We have studied the asymptotical behaviour of the solutions when t — > 00, finding that there is a simple asymptot- 
ically stable anisotropic solution for all cases. It can be easily shown (see [6|) that this solution remains anisotropic 
for large t. 
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